Spectral collapse in ensembles of meta-molecules 
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We report on a new collective phenomenon in metamaterials: spectral line collapse with increasing 
number of the unit cell resonators (meta-molecules). Resembling the behaviour of exotic states of 
matter, such as Bose-Einstein condensates of excitons and magnons, this new effect is linked to 
the suppression of radiation losses in periodic arrays. We demonstrate experimentally spectral line 
collapse at microwave, terahertz and optical frequencies. It emerges as a universal and truly scalable 
effect underpinned by classical electromagnetic interactions between the excited meta-molecules. 



The burgeoning field of metamaterials provides unique 
opportunities to engineer the electromagnetic properties 
of artificial media and achieve exotic functionalities, such 
as negative refraction [l| and cloaking 0. Similarly to 
natural crystals, which are created by arranging indi- 
vidual atoms and molecules in a regular grid, periodic 
ensembles of subwavelength electromagnetic resonators 
present an effective medium to an incident with proper- 
ties not available in natural materials. Here, we study 
the dependence of the metamaterial properties on the 
number of meta-molecules in the microwave, THz and 
optical domain, and demonstrate a new collective phe- 
nomenon in metamaterials: in contrast to solid state 
crystals, where bulk arrangements result in broadening of 
the individual clement spectral line, leading eventually to 
the formation of absorption bands, regular ensembles of 
meta-molecules can exhibit the opposite effect, i.e. spec- 
tral line collapse. 
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FIG. 1: (Colour online) Resonant absorption lines of atoms 
form a broader absorption band in solids. In coherent meta- 
materials we see the opposite trend: resonant lines of large 
arrays are much narrower than lines of individual meta- 
molecules. 

The reported phenomenon is characteristic to a novel 
class of artificial media, which we call " coherent" meta- 



materials [3j and are characterized by very strong in- 
teractions between the elcctromagnctically excited meta- 
molecules that provide for a low rate of energy loss due to 
scattering and lead to a high-quality resonant response. 
An example of a coherent metamaterial is an array of 
ASRs, where the meta-molecular excitation corresponds 
to an oscillating magnetic dipole perpendicular to the 
plane of the array that does not interact directly with 
the magnetic field of the incident wave, thus creating a 
nearly thermodynamically isolated ensemble of strongly 
interacting coherent "molecules" with interesting physi- 
cal properties. To illustrate this behavior we present a 
comparison with an "incoherent" metamaterial: a two- 
dimensional array formed by pairs of concentric conduct- 
ing rings that also supports a high-quality resonant re- 
sponse. In this case, however, the response of the array is 
a sum of the individual meta-molecule responses, rather 
than a collective property. 

The coherent microwave metamaterial was manufac- 
tured as a regular planar array of asymmetric split rings 
(ASR) etched from a 35 ^m thick copper layer on a 
1.6 mm thick FR4 substrate. The diameter of the ASR 
was 6 mm with a line width of 0.4 mm and was split 
in two segments corresponding to 140° and 160° arcs. 
The unit cell of 7.5 x 7.5 mm 2 rendered the arrays non- 
diffracting at normal incidence for frequencies of up to 
40 GHz. In the incoherent metamaterial ASRs were re- 
placed with pairs of concentric rings. The inner and outer 
rings had the diameters of correspondingly 4.50 mm and 
5.45 mm, and were both 0.2 mm wide. The transmission 
measurements were performed in a microwave anechoic 
chamber at normal incidence using broadband linearly 
polarized horn antennas equipped with collimating lenses 
and a vector network analyzer. 

Transmission spectra of large coherent and incoherent 
arrays show similar resonant features in the form of a 
broad stop-band split by a narrower Fano-likc transmis- 
sion peak (see Fig. 2). For the array of ASRs the trans- 
mission resonance occurs at around 11 GHz and is as- 
sociated with the excitation of anti-symmetric currents 
oscillating in the opposite arcs of each split ring (see in- 
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FIG. 2: (Colour online) Transmission spectra feature trapped- 
mode resonances of the coherent (asymmetric split rings, 
solid red curve, v co h) and incoherent (concentric rings, solid 
blue curve, fine) metamaterial arrays of 88 meta-molecules. 
Dashed curves indicate transmission stop bands of the meta- 
materials in the absence of trapped-mode excitations. The 
top right inset shows a sample of 88 meta-molecules (unit 
cells) of the coherent array exposed through a metal mask 
(represented by the darker area). Other insets show the 
trapped current modes excited in the unit cells of the co- 
herent (left) and incoherent (right) metamaterials and the 
corresponding induced magnetic dipole moments m. 

set to Fig. 2). Such a current mode, known as a trapped 
or closed mode, yields magnetic dipole moments oscil- 
lating synchronously (coherently) in all meta-molecules 
along the direction normal to the plane of the array. The 
induced magnetic dipoles interact strongly with one an- 
other, while their interaction with the perpendicularly 
oriented magnetic field of the incident electromagnetic 
wave is forbidden For the double-ring metamaterial 
a similar narrow resonance at 13 GHz corresponds to os- 
cillations of oppositely directed currents in the inner and 
outer rings, as shown in the inset to Fig. 2. Although such 
current configurations give rise to magnetic moments, 
the latter cancel one another and the total magnetic mo- 
ment of the meta-molecule and thus interactions between 
meta-molecules are negligible [H[ . 

Figure 3 shows the dependencies of the transmission 
resonance quality factors, Q = v/Ais, on the total num- 
ber of meta-molecules, N, that form the arrays. In the 
experiment, the number of meta-molecules exposed to 
the incident wave was controlled by placing metal masks 
of different sizes on a large metamaterial array, as illus- 
trated in the inset to Fig. 2. These masks screened the 
peripheral parts of the array, leaving the central part 
of the array exposed to electromagnetic radiation. The 
shape of the masks was close to circular with a step-like 
profile of the opening, which ensured that every unit cell 
of the array was either fully screened or exposed. These 
measurements were performed with 22 different masks 
exposing arrays with a total number of unit cells in the 



range from 32 to about 700. Experiments with larger 
arrays were not practical and unnecessary as already for 
N ~ 700 the quality factor saturated. Experiments with 
smaller arrays were not feasible as for N < 32 diffraction 
takes its toll on the accuracy of the data. In the mea- 
sured range of N the experimental data clearly indicate 
that the Q- factor of the coherent metamaterial strongly 
depends on the size of the illuminated area, i.e., the to- 
tal number of meta-molecules engaged in the interaction 
with the incident wave. Indeed, the Q-factor, which mea- 
sures about 10.5 for the smallest array (N = 32), can be 
seen to increase by almost 70% reaching 17.5 for the full- 
sized array with N = 688. In contrast, the Q-factor of the 
incoherent metamaterial appears to be practically inde- 
pendent of the number of exposed unit cells and remains 
at around Q = 9. 
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FIG. 3: (Colour online) Graphs show Q-factors as a function 
of the total number of meta-molecules. Experimentally mea- 
sured data are represented by points for both coherent (di- 
amonds) and incoherent (circles) metamaterials. Solid lines 
present theoretical fits to the data. 

To illustrate further the dependence of the resonant 
response on the size of the coherent metamaterial, we 
studied the spatial distribution of the magnetic field in a 
single split-ring meta-molecule and a meta-molecule that 
is a part of a regular array containing TV = 400 unit 
cells. For this we employed a THz near-field imaging 
technique with sub-wavelength resolution described in 
Q. The technique enabled accurate mapping of orthog- 
onal components of the electrical field E in the plane of 
the unit cell (as illustrated in the inset to Fig. 4a), which 
was used to calculate the component of the magnetic field 
H normal to the plane of the structure via the Maxwell 
equation V x E = —dB/dt. 

In our experiment we observed that at resonance con- 
ditions the spatial distribution of the magnetic field was 
radically different for a split-ring placed in the array and 
an isolated split-ring. For a ring in the array the mag- 
netic field penetrates the unit cell in the same direction 
everywhere within the area enclosed by the ring: the 
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FIG. 4: (Colour online) Panel (a) shows the spectral depen- 
dence of the squared total magnetic flux amplitude through 
an ASR for the cases of one isolated ring (empty squares, 
weak magnetic response) and a ring located in the middle of a 
20 x 20 array (filled squares, strong magnetic response). Pan- 
els (b) and (c) present snap-shot maps of the instantaneous 
component of the magnetic field normal to the plane of the 
ring (colour coded) and the in-plane component of the electric 
field (arrows) at resonance frequencies v\ and V2 correspond- 
ingly (see also supplementary videos for temporal evolution 
of the fields). In order to show both fields at their maxima on 
the same map the magnetic and electric fields are presented 
with phase shift of tt/4. 

total magnetic moment of the meta-molecule is at max- 
imum indicating the excitation of strong anti-symmetric 
(i.e. ring) currents in opposite sections of the ring (see 
Fig. 4b). while the spectral dependence of the squared 
amplitude of the magnetic flux reveals a narrow sharp 
peak centered at v\ — 0.165 THz (see Fig. 4a). For an 
isolated single meta-molecule the net magnetic flux is 
also non-zero but it is much lower than for a ring in the 
array: the magnetic field is oppositely directed in adja- 
cent sections of the ring (see Fig. 4c), the total magnetic 
flux is small and the anti-symmetric current component is 
weak. As a result, the resonant feature at v% = 0.142 THz 
is much less pronounced and its Q-factor is significantly 
lower than for rings in the array indicating considerable 
damping (see Fig. 4a). 

We attribute the observed size-dependent resonant re- 
sponse of the coherent metamaterial to the existence of 
strong interactions between the meta-molecules mediated 
by magneto- inductive surface waves 0. It follows from 



symmetry arguments that magnetic dipoles coherently 
oscillating perpendicular to the plane of an infinite (or 
very large) regular array can re-radiate only in its plane, 
which ensures very low scattering losses and therefore 
a high Q-factor of the system. Electromagnetic energy 
in the coherent array is trapped in the form of standing 
magneto-inductive surface waves, which may only scatter 
on the edges of the array. The scattering losses grow and 
the Q-factor diminishes as the array becomes smaller, 
and in the limiting case of an isolated meta-molecule the 
resonance is extremely weak due to intense energy losses 
through magnetic dipolc radiation (see also supplemen- 
tary materials). 

The dependence of the Q-factor on the size of the ar- 
ray can be qualitatively understood from the definition 
of the Q-factor, Q = 2ir ■ E/AE, where E is the energy 
stored in the array and AE is the energy dissipated by 
the array per cycle. Here E is simply proportional to 
N, while AE = AE m + AE S represents contributions 
from two loss mechanisms. AE m = a m ■ N includes ma- 
terial losses, i.e., Ohmic losses in the metal and dissipa- 
tion losses in the dielectric substrate, as well as radia- 
tion losses due to an electric dipole induced in each unit 
cell, which is controlled by the degree of asymmetry of 
the split ring. Losses associated with scattering of the 
magneto- inductive waves on the edges of the array, AE S , 
are proportional to the length of its perimeter and thus 
AE S = a s -y/~N. Thus, the phenomcnological dependence 
of the quality factor on the number of unit cells can be 
presented in the form Q cx Nj (a m -N + <j s \/~N), which fits 
the experimental data very well (solid red line in Fig. 3). 
In comparison, in the case of the incoherent array of con- 
centric pairs of rings, the total induced magnetic moment 
perpendicular to the plane of the array is small, hence 
magnetic dipole interactions and scattering of magneto- 
inductive waves are negligible (a m ■ N 3> a s vN). The 
second term in the denominator can be disregarded mak- 
ing the Q-factor independent of the size of the array, 
which is in complete agreement with our experimental 
observations (solid blue line in Fig. 3). 

The dependence of the resonant response on the ar- 
ray's size results from truly classical interactions between 
the excited states of the meta-molecules and therefore 
shall be universal and truly scalable with wavelength. 
Since the characteristic size of the array at which sub- 
stantial broadening of the resonance is seen is given by 
N = (er s /<7 m ) 2 , with increase of Joule losses the size de- 
pendent effect will be seen in smaller arrays, i.e., in an 
array with higher Joule losses the coherent state will be 
formed by a smaller number of meta-molecules. 

This is exactly what we observed in our experiments 
in the optical part of the spectrum, where Joule losses 
are dominant. Here we studied transmission of square 
asymmetric split ring slit arrays of different sizes scaled 
down 15,000 times relative to the microwave sample. The 
ring slits (complimentary to split-rings) were cut from 
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a 55 nm thick gold film supported by a 500 /xm silica 
substrate using focused ion beam milling. The transmis- 
sion measurements were performed using a microscope- 
based spectrophotometer by CRAIC Technologies. Be- 
cause of the negative configuration the size of the illu- 
minated area in transmission was naturally controlled by 
the apertures of the arrays, while unstructured square 
windows of the same size were used as reference. We 
studied arrays of five different sizes containing from 16 
to 144 meta-molecules (see Fig. 5a - 5b). Such pho- 
tonic mctamaterials exhibited a trapped-modc resonance 
at around 1000 nm, which, in full accordance with Babi- 
net's principle, was seen not as a transmission peak (as 
for a positive array, see Fig. 2), but as a dip in transmis- 
sion for the largest array of 12 x 12 unit cells (see inset 
to Fig. 5c). For smaller arrays the resonance gradually 
became broader and shallower and, as shown in Fig. 5c, 
completely disappeared for N = 16, which suggests that 
the most dramatic change occurs when the number of 
unit cells drops below 49. 
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FIG. 5: (Colour online) Panel (a) shows SEM images of coher- 
ent metamaterial arrays of different sizes. Panel (b) presents 
a schematic of the metamaterial unit cell. Panel (c) shows 
transmission spectra of the metamaterial arrays in the vicin- 
ity of the trapped-mode resonance. The inset presents the 
transmission spectrum of the 12 x 12 array over a much wider 
wavelength range; a dashed box indicates the spectral domain 
that is covered by the main plot, while an arrow marks the 
position of the trapped-mode resonance. 



The results presented above illustrate that larger ar- 
rays of coherent mctamaterials exhibit narrower reso- 
nances. They also show that in samples with smaller 
losses more meta-molecules are engaged in forming the 
coherent state: the Q-factor starts reducing dramatically 
below about 200 meta-molecules in the microwave ar- 
ray, while in the lossy optical case we saw a significant 
reduction of Q-factor below 49 meta-molecules. We ar- 
gue that the number of molecules forming the coherent 
response will increase as Joule and dielectric losses are 
compensated. This is relevant to the recent observation 
of loss compensation in a coherent photonic metamate- 
rial [|| , which paves the way to an intriguing opportunity 
of creating a lasing spaser, a coherent source of optical 
radiation fuelled by coherent plasmonic oscillations in in- 
dividual metamaterial resonators of a coherent array. In 
the latter case the gain substrate supporting the rings 
would be the source of energy Q . 

In conclusion, we would like to point out that the 
spectral collapse of the metamaterial resonance with in- 
creasing number of meta-molecules has intriguing phe- 
nomenological similarity to the narrowing of lumines- 
cence/scattering spectra upon formation of the Bosc- 
Einstein condensates (BEC) of exciton and magnon in 
solids [ljj [ll], G3 ■ The latter is observed when concen- 
tration of the bosonic quasi-particles in equilibrium ex- 
ceeds a certain critical value and most of the particles are 
statistically forced to occupy the ground level forming a 
stable coherent eigenstate. In the case of the excitons, for 
example, the narrow emission line of the condensate cor- 
responds to the binding energy of the quasi-particles and 
results from their recombination in the condensate. The 
analogy between these two phenomena becomes appar- 
ent since the ground state of the metamaterial system 
is achieved through de-excitation of the low scattering 
(i.e. trapped) electromagnetic mode of interacting meta- 
molecules. Moreover, as in the BEC, the collapse of the 
resonant line is most rapid when the number of meta- 
molecules reaches a certain critical value (N c as 100 in 
the case of the microwave metamaterial). 
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